Exchange induced charge inhomogeneities in rippled neutral graphene 
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A new mechanism that induces charge density variations in corrugated graphene is proposed. 
Here it is shown how the interplay between lattice deformations and exchange interactions can 
induce charge separation, i.e., puddles of electrons and holes, for realistic deformation values of 
the graphene sheet. The induced charge density lies in the range of 10 11 - 10 12 cm , which is 
compatible with recent measurements. 

PACS numbers: 73.21.-b,73.20.Hb,73.22-f 



o 
o 



> 

oo 
in 
m 

(N 
O 

o 



X 



Introduction. The ability to isolate and manipulate 
graphene -single layers of carbon atoms packed in a hon- 
eycomb lattice- has opened and boosted the experimental 
study of the properties of two dimensional (2D) mass- 
less Dirac fermionsp], [H, Q. The existence of a strictly 
2D crystal structure is puzzling in itself, as, according to 
theory and experiments, perfect 2D crystals can not exist 
in the free state 0, H, 0]. 

Recently, transmission electron microscopy 
experiments [7| have resolved that suspended graphene 
sheets are not perfectly flat but exhibit microscopic 
roughness or ripples such that the surface normal varies 
by several degrees and the out of plane deformation can 
reach lnm. This deformation corresponds to a rather 
large strain, 1% , but is sustainable without plastic 
deformations or generation of defectsQ- Also, layers 
placed on SiC>2 seems to follow the corrugation of the 
substrate [1, Q. The height and width of these ripples 
are consistent with models which allow the carbon ions 
to form different types of bonds [Toj. 

Recent single electron transistor based experimentsflU 
have evidenced the existence of electron and hole pud- 
dles of densities ~ 10 10 — 10 11 cm -2 in the vicinity of the 
neutrality point. The existence of these puddles could be 
simply related to the presence of a disorder potential in- 
duced by the substrate (T3, [l3j . An alternative to this ex- 
planation bears on the existence of ripples which leads to 
a modulation of the hopping amplitudes between carbon 
atoms. This modulation affects the electronic structure 
in two-fold manner. Firstly, the modulation induces an 
effective magnetic field which changes locally the density 
of states, but does not separate charge In this respect 
it has been argued in Ref. [IE] that a one-dimensional 
deformation of the graphene sheet will form zero energy 
Landau levels corresponding to an effective magnetic field 
of tens of Teslas. This will increase the compressibility of 
the system and eventually would induce electronic phase 
separation. Secondly, second-neighbor hopping changes 
can induce a potential on the carriers which does sepa- 
rate charge 1J] . An additional source of change in the 
local density of states has been suggested to come from 
local Fermi velocity changes induced by the curvature 
associated to the ripples [16}. This effect induces charge 




FIG. 1: (Color online) Schematic view of graphene with pe- 
riodic ripples. L denotes the lateral extent of the ripples and 
h the out of plane displacements of the Carbon atoms. 



inhomogeneities in doped graphene, but, in the presence 
of electron-hole symmetry, can not explain the existence 
of electron-hole puddles in undoped graphene. 

In this paper we propose that because of the depen- 
dence of the exchange energy on the density of carbon 
atoms, the strain modulation produced by the ripples 
induces a charge modulation in undoped graphene. In 
order to emphasize and highlight the effect of the ex- 
change on the charge inhomogeneity we have disregarded 
in this work the dependence of the hopping parameter 
on next neighbor distance. This is a justified assump- 
tion as in graphene a variation of the hopping amplitude 
can be described by means of a fictitious vector potential 
which, in general, does not induce significant changes in 
the density of statesJTij]. In other words, the degener- 
acy of the Landau levels and their contribution to the 
compressibility [TH will be strongly suppressed if a ran- 
domly curved graphene sheet fK, [l0| is considered. 

In the following we first point out how the existence 
of ripples in the graphene surface creates a modulation 
of the distance between first neighbor carbon atoms. 
Then, within the Hartree-Fock approximation, we obtain 
the total energy of the system as a function of the dis- 
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tance between first neighbors atoms, d, and of the extra 
charge per carbon atom, <5. Finally, we show that a long- 
wavelength modulation of the lattice parameter in the 
graphene sheet can induce a charge density modulation 
in undoped graphene. 

Distance between carbon atoms. To simulate the ex- 
perimental electron diffraction patterns Q, it is needed 
to assume a ratio L/h w 10 between the lateral size of 
the ripples, L, and the out of plane displacement, h (see 
Fig. A prudent estimate for the typical lateral size of 
the crumpling has been estimated to be between 2 and 
20 nm 7]. In our calculation wc consider the following 
out of phase modulation, 



z(Ri) = h [sin {GR X . % ) + sin (G R Vti ) 



(1) 



where G = 2ir/L and the vectors Ri are the position 
of the carbon atoms in a perfectly flat graphene sheet. 
The average distance between an atom and its three first 
neighbors depends on the position of the considered atom 
and is given by 



d(Ri) ~ d 1 



h 2 G 2 
8 



[cos (2GR x ,i) + cos (2 G R V}1) 



(2) 

being d = do(l + h 2 G 2 /A) the average distance be- 
tween carbon atoms in the presence of the ripple and 
do the equilibrium distance between the C atoms in flat 
graphene. Expression @ indicates that a modulation in 
the out of plane position of the carbon atoms implies a 
modulation in the average distance between them. 

Electronic structure and kinetic energy. The electronic 
active states of graphene are the bands formed by the 
carbon p z orbitals which are ordered in a honeycomb lat- 
tice (triangular lattice with two atoms per unit cell). The 
band structure is accurately described by means of a first 
neighbors tight-binding Hamiltonian with a unique hop- 
ping parameter t. Because the bipartite character of the 
honeycomb lattice, the hamiltonian in reciprocal space is 
represented by a 2x2 matrix. In undoped graphene there 
is one electron per carbon atom and the conduction and 
valence bands touch at two non-equivalent points of the 
Brillouin zone that are called Dirac points. 

For wavevectors near the Dirac points the electronic 
structure can be described by a Dirac Hamiltonian of 
the form 



Hkin = ftv F k • a 



(3) 



where a are the Pauli matrices, k is the electron mo- 
mentum measured with respect the Dirac point and the 
Fermi velocity only depends on the hopping energy t and 
the next neighbors distance do, 



hv F = - t d 



(4) 



The eigenvalues and eigenvectors of Hamiltonian Eq[3]are 



e ikr / 1 

e±(k) = ±hv F \k\ and V±(k) = —j= I ±pU 



, (5) 



where 9^ = arctanfc^/fcj,. 

The use of the continuum Dirac Hamiltonian for de- 
scribing the properties of graphene requires the introduc- 
tion of a maximum value of momenta, k c , which is chosen 
to keep the number of states in the Brillouin zone fixed, 
i.e., gditk 2 = (2tt) 2 /Sq. Here = 2 is the Dirac points 
degeneracy and So = 3\/3d 2 j /2 is the unit cell area. 

When doping graphene with electrons or holes the 
extra carriers form a Fermi sea with Fermi wavevector 
kp = y/Aii\n\/ {gd9sj, where n is the 2D density of added 
charge and g s = 2 is the spin degeneracy. The kinetic 
energy of the system per carbon atom is given by 
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where 8 is the extra charge per carbon atom with respect 
the intrinsic situation. Note that due to the linear dis- 
persion of the bands, when we discard the variation of 
t on distance, the kinetic energy per carbon atom does 
not depend on the distance between atoms. From the 
expression of the kinetic energy, Eq. ([6]) , there is always a 
kinetic energy cost associated to modulating the charge 
in undoped graphene. 

Exchange energy. The exchange contribution to the 
total energy per carbon atom has the form[l7l [l8l. [l9j 



E 



ex = -^So Kg)^ lS '(k,k+qX(k)n s ,(k+q), 

s,s' ,k,q 

(7) 

where s and s' is the band index (±1), v(q) = 2ire 2 /eq 
the 2D Fourier transform of the Coulomb interaction, e is 
the dielectric constant of the system, n s (k) is the Fermi 
occupation function of the state {s, k} and F s>s i (k, k+q) 
is the square of the overlap between the wavefunctions, 
tp s (k) and Vv(k + q) 
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F S)S ,(k,k + q) 
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with 9 the angle between the wavevectors k and k + q. 
For Coulomb interaction, the factor ((8]), makes the ex- 
change interaction to be larger between states in the same 
band. In the expression Eq. ([7]) we have neglected the 
exchange energy between electrons belonging to different 
Dirac cones. This is appropriated in the long wavelength 
approximation because the difference in momentum be- 
tween states coming from different Dirac points is very 
large. 

Following the notation of Ref. [13 the exchange energy 
per carbon atom can be written as 



E EX (d, S) 
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where the functions R n (a) are defined in Ref. In the 
limit of small extra charge one can approximate 



(10) 



In Eq. [5] the first term corresponds to the exchange en- 
ergy of the full occupied valence band, the second term 
the exchange energy of the extra electrons or holes, and 
the last term is the variation of the exchange energy be- 
cause the interaction between the extra carriers and the 
valence band electrons. Note that the last term in Eq. [9] 
changes sign with the electron/hole character of the extra 
carriers. 

Coupling between lattice deformation and charge. Ex- 
pressions Eq. ([6]) and Eq. ([9]) have been obtained for a 
uniform system. For long wavelength modulations of d 
and S, the total energy per carbon atom can be written 
as 



Eh 



dr[E kin (6(r)) + E E x{d(r),6(r))} . (11) 



This equation couples a modulation of the first neighbors 
distance between carbon atoms with a modulation of the 
electron density. 

In undoped graphene a modulation d(R) of the form 
described in Eq[2] will produce, in lineal regime, a mod- 
ulation of the electric charge of the form, 



S(R) = Si [cos (2 G Rxt) + cos (2 G R y 



(12) 



where the amplitude of the charge modulation, Si, is ob- 
tained by minimizing the total energy Eq llll To lowest 
order, the charge amplitude Si only depends on the ratio 
h/L and on a dimensionless coupling constant defined as 
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e 2 Ao 



(13) 



The constant g indicates the ratio between the Coulomb 
and the kinetic energy of the electron system. 

For values larger than a critical coupling constant 
g c ~ 2.18, the total energy does not present a minimum 
as function of St. This indicates that the system is un- 
stable against more correlated phases and therefore the 
lower order calculation fails. This value of the coupling 
constant is considerably smaller than the critical value 
of g (g ~ 5.3) needed for the occurrence of instabilities 
towards highly correlated states in flat graphene [lTf. For 
values of g < 2.18 the total energy presents a minimum 
and the electronic system reacts to the lattice deforma- 
tion by modulating the charge according to Eq. 1121 In 
Fig|2] we plot the value of Si that minimize the total en- 
ergy as function of the coupling constant for different 
values of the ratio h/L. 

The distance between carbon atoms in graphene and 
the hopping amplitude are quantities well established; 
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FIG. 2: (Color online) Amplitude of the charge modulation 
per Carbon atom Si as function of the coupling constant, for 
different values of the ratio h/L. 



t = 2.5eV and do=lA2A. Therefore a realistic estima- 
tion of the coupling constant, g, only requires the knowl- 
edge of the value of the dielectric constant. This should 
take account both the effect of the screening current in 
the substrate as well as the weak intrinsic screening in 
graphene [13] • We estimate that depending whether the 
graphene sheet is freely suspended or it is placed on top 
of an insulator, the dielectric constant would take values 
in the range 3 < e < 4. For these values of e, the coupling 
constant takes values in the range 0.66 < g < 0.88. For 
these values of g, we can see in FigJ5]that the use of the 
lower order coupling between h and Si is appropriated. 

For an intermediate coupling constant, g — 0.75, we 

obtain Si ~ 9 (jj^s) which, for values of the ratio h/L 

in the range 0.05 to 0.1, implies values of Si in the range 
5.6x 10 -5 to 0.9x 10 -4 . This corresponds to density mod- 
ulations in the range 2.25xlO n cm~ 2 to 3.6xl0 12 cm~ 2 . 
From this estimate we conclude that a modulation of the 
out of plane position of the carbon atoms of amplitude 
h ~ 1 — 2 nm in a lateral size L ~ 10 — 20nm induces 
a modulation in the charge density of order 10 11 cm" 2 . 
This magnitude of the charge modulation agrees with the 
density of charge in the electron-hole puddles observed by 
single electron transistor based experiments [ll| . Interest- 
ingly, in suspended graphene the contribution of the sub- 
strate to the dielectric constant is practically suppressed. 
In this case the coupling constant becomes larger than 
in the case of graphene placed on a dielectric. There- 
fore we expect a higher density modulation in suspended 
graphene than the observed in graphene on Si02- 



It is also pertinent to estimate the importance of the 
Hartree repulsion on the values of the charge density 
modulation. The Hartree energy per carbon atom takes 
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the form 

* - 

G' 

8tt 2 3\/3 ^edo 

In the above expression n(G') is the G' component 
Fourier transform of the charge. For the values of dielec- 
tric constant and L we obtain that the Hartree energy is 
much smaller than the kinetic and exchange energy and 
the values obtained for the charge density modulation are 
practically unaffected by the Hartree repulsion. 

On the top of the effect discussed in this work one 
should also consider the influence of the hopping depen- 
dence on lattice strain discussed in Refs. [lj, [l5l Realis- 
tic estimates of the influence of second-neighbor hopping 
variations indicate that these are typically smaller that 
the one discussed here. In any case, both effects are 
compatible. Regarding the formation of Landau levels 
and their influence on the compressibility, we have veri- 
fied numerically 23] that that density of states close the 
neutrality point is not seriously affected by the effective 



magnetic fields when crumpled graphene with uncorre- 
cted ripples is taken into consideration. This contrasts 
with the strong influence that periodic ripples exert on 
the density of states near the Dirac point [Fa]. 

Summary. We notice that exchange interaction be- 
tween carriers produces a coupling between a modulation 
of the distance between first neighbors atoms in graphene 
and a charge density modulation. This mechanism con- 
nects the presence of ripples in undoped graphene with 
the existence of electron hole puddles of density up 
to 10 12 cm -2 . These densities are of the same order 
than the observed recently in single electron transistor 
measurements [ill- In suspended graphene the absence 
of substrate makes the exchange energy stronger and we 
predict that this increase will produce a higher modula- 
tion of the charge density. 
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